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In this paper, we introduce a Ginzburg-Landau (GL) theory for the extended-s and d-wave super- 
conductors (SC) in granular systems that is defined on a lattice. In contrast to the ordinary Abelian 
Higgs model (AHM) that is a GL theory for the s-wave SC, Cooper-pair field (Higgs field) is put on 
links of the lattice in the present model. By means of Monte-Carlo (MC) simulations, we study phase 
structure, gauge-boson mass (the inverse magnetic penetration depth) and density of instantons. In 
the ordinary noncomapct U(l) AHM, there exists a second-order phase transition from the normal 
to SC states and the gauge-boson mass develops continuously from the phase transition point. In 
the present gauge system with link Higgs field, on the other hand, phase transition to the SC state 
is of first order at moderate coupling constants. The gauge-boson mass changes from vanishing to 
finite values discontinuously at the phase transition points. 



Ginzburg-Landau (GL) theory plays a very important 
role for study on superconducting (SC) phase transition. 
By taking into account the effect of the electromagnetic 
interactions, one can show that it has a form of a noncom- 
pact U(l) lattice-gauge- Higgs model in which the order- 
parameter boson field sits on lattice sites [1]. In the last 
few decades, unconventional superconductors, whose or- 
der parameter of Bose condensation is not the usual s- 
wave, have been discovered [2]. In this paper we shall 
extend the GL theory to that for unconventional SC like 
the d x 2_ y 2-w&ve SC in which the order-parameter field, 
the Cooper-pair wave function, changes its sign under a 
ir/2 rotation of the real-space coordinates. Therefore in 
order to describe such an unconventional SC, the order- 
parameter field, Higgs boson field, must be put on lat- 
tice links instead of lattice sites. On-site amplitude of 
the Cooper pair is vanishingly small because of, e.g., the 
strong on-site Coulomb repulsion. 

We shall use the path-integral formalism, and define 
the model on a three-dimensional (3D) cubic lattice of 
system size L? with the periodic boundary condition. Be- 
fore going into details of the extended model, let us first 
consider the ordinary noncompact Abelian Higgs model 
(AHM) as a reference system, which is defined by the 
following action, 

A A hm = x [X] CuF ?A x ) + X K <t>x+i u xj<l>l > (!) 

pi link 

where i^(ar) = A Xyi - A x+j4 + A x+itj - A xJ = 
1,2,3) and a gauge field A x ^ on the link (x, i) is re- 
lated to the electromagnetic vector potential A? m as 

4m = & m ■ and U x,j = z lA *- ] - 4>x is the Higgs 
field corresponding to the s-wave Cooper pair and in the 
London limit 4> x = e t(Pm (tp x £ [—w,tt]). c^ 1 is the elec- 
tric charge of the Cooper pair and k is a parameter corre- 
sponding to the superfluid density and a decreasing func- 
tion of the temperature (T). Aahm is nothing but the 



3D XY model coupled with the noncompact U(l) gauge 
field that describes the electromagnetic interactions. 

We studied the phase structure of the model Aahm 
by the Monte-Carlo (MC) simulations calculating the in- 
ternal energy E = — (Aahm)/^ 3 and the specific heat 
C = ((Aahm — (Aahm)) 2 ) / L 3 ■ We found that there ex- 
ists a second-order phase transition line [3] emanating 
from the 3D XY critical point at (k, c u ) = (0.46, oo) (ob- 
tained in the system of size L = 24). In Fig. 1(a), we show 
the specific heat as a function of k with c u = 1. The spe- 
cific heat exhibits a typical behavior of the second-order 
phase transition as the system size L 3 is increased from 
8 3 to 24 3 . This result is in sharp contrast to the compact 
AHM in the London limit in which no phase transitions 
occur and only the confinement phase exists [4]. 

We also measured the gauge-invariant gauge-boson 
mass Mq which is defined through the correlation func- 
tion of the operator sin(-Fy(;r)). More precisely let us 
define the operator 0{x) as, 

0(x) = X eij8inFij(x), eia = -£21 = 1, (2) 

i,i=l,2 

and its Fourier transformed operator in the 1 — 2 plane, 
6(2:3) = 0(x)e lpiXl+lp2X2 . (3) 

One can expect that the correlation function of 0(^3) 
behaves as 

(d(x 3 )0(x 3 + *)) cx e-Vpf+Pl+^c *. (4) 

In Fig. 1(b), we show the result. We define the gauge- 
boson mass Mq from numerical results as Mq = 
sign (A 2 — p 2 )i/A 2 — p 2 , where A is the inverse correlation 
length of the Fourier transformed operator 0(^3) with fi- 
nite momentum p [5]. The negative values for k < 0.52 
in Fig. 1(b) come from the above definition of the mass 
Mq and it is considered as a finite-size effect [5] . 
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From Fig. 1(b), we conclude that Mq is vanishing for 
k < 0.52 and develops continuously as n increases. This 
behavior obviously is consistent with the specific heat 
measurement in Fig. 1(a) and indicates the existence of 
a second-order phase transition from the normal to the 
Higgs-SC phases, though the value of the critical coupling 
k c obtained from Mq (k c = 0.52 with L = 16) is slightly 
different from that obtained by C (k c — 0.54 with L = 
16). Similar phenomenon has been observed also in the 
previous studies on the U(l) gauge field coupled with the 
CP 1 fields [6]. 
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FIG. 1. Results of the MC simulations for the AHM with 
c-u = 1. (a) Specific heat C as a function of k. It in- 
dicates the existence of a second-order phase transition at 
jc c = 0.538 (L = 24). (b) Gauge-boson mass Mq obtained 
from the correlation function of sin(_FV, (2;)). The critical cou- 
pling is estimated as k c = 0.52 (L — 16). 

We applied the finite-size scaling to C as C(k,L) = 
L a I v §(1?- 1 v e) , e = (k-KkjJ/kqo, where is the critical 
coupling at L — > oo and 4>(x) is the scaling function. We 
found that — 0.534, v — 0.83 and a = 0.15. 

Let us introduce the extended mode that is defined by 
the following action, 

Ag ^ = \zZ [ c ^ f U x ) + c - yi + c ^(uvuv + VUVU) 



pi 



+d m (UUVV + 3 cyclic permutations) 



(5) 



where V x j is the spin-singlet Cooper-pair field on link 
that is related to electron operator ip xcr (cr =f, J.) as 



V x<j oc {ip x fip x+ji - ip xl ip x+j] ). 



(G) 



From Eq.(6), it is obvious that V x> j can be regarded as 
a three-component vector field. Here it is interesting to 
notice that a GL theory for the spin-triplet p-wave super- 
conductivity in ferromagnetic ZrZn2 was proposed and it 
employs a SC order parameter similar to V x j [7]. Gradi- 
ent terms of the GL theory for ZrZn2 have a similar form 
to ^GLin Eq.(5). 

Hereafter we shall consider the London limit of V x j 
and set V x j — e l8wj (6 X j G [— n, it]). Each term in Aql 
is depicted in Fig. 2 where c u , c v etc are coupling con- 
stants, and Aql is constructed to be invariant under the 
following noncompact local gauge transformation, 
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(7) 



From Eq.(7), it is obvious that V x j can be regarded as 
another gauge field dual to the electromagnetic gauge 
field [8]. We consider terms as local as possible for Aql, 
and the partition function Z is given as 



[DA] / [DO] e 



(8) 
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FIG. 2. Action A GIj of the GL theory (5). 

Compact U(l) version of the above Aql, in which 
J2„i Ffj{ x ) is replaced by : f/ 4 , has been studied in 



the previous paper [9]. In the present paper, we study 
the noncompact U(l) gauge theory as a GL theory for 
the unconventional SC in which the gauge field U x j de- 
scribes the electromagnetic field. 

There is credible evidence that the SC phase transition 
in the high-T c cuprates is of second-order and further- 
more it is in the 3D XY model universality class [10]. On 
the other hand as explained above, the Cooper-pair field 
must sit on lattice links instead of sites in order to de- 
scribe the d-wave SC state. One of the simplest GL the- 
ory for the d-wave SC is Aql given in (5) . In the present 
paper, we shall study the phase structure and physical 
properties of Aql by means of the MC simulations and 
compare the results with those of the noncompact AHM 
and the XY model. 

Let us first notice that for vanishing c m — d m = 0, 
the present system reduces to two independent decoupled 
gauge models, a noncompact U(l) gauge model of A x j 
and a compact U(l) gauge model of V x j. Then it is 
obvious that there exist no phase transitions in that case. 
By giving finite values for c m and/or d m , we study phase 
structure in the c u — c v plane. 

We first study the case c m = 1, d m — 0. It is in- 
structive to consider the large-c u limit in which config- 
urations of the noncompact gauge field are restricted as 



.4 



x,3 



'■P. 



x+j 



tp x . Then the c m -terms in Aql become as 



C.C. 



(9) 



We shall call the above term double Higgs coupling. As 
we explained above, the usual Higgs coupling of the com- 
pact gauge field V x j, J2x e.~ lip=c V X je~ l ' f '* +: > does not 



induce any phase transition. On the other hand, the 
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doubly-charged Higgs coupling, e~ lv * (V Xl i) 2 e~' l<Pm + i , in- 
duces a phase transition from the confinement to Higgs 
phases [11]. Then it is interesting to study the extended 
model also from the viewpoint of the Higgs coupling and 
see if the Higgs phase transition occurs as the double 
Higgs coupling (9) is increased. 

We studied the phase structure in the c u — c v plane 
by calculating E and C and found that there is a phase 
transition line. Typical behavior of E near phase tran- 
sition points is shown in Fig. 3 (a), which indicates that 
the transition is of first-order. To see physical meaning 
of the phase transition, we measured the instanton den- 
sities of the gauge fields U x j and V x j. We follow the 
definition of the instanton densities pu and py given in 
Refs. [12,9]. The result is shown in Fig.3(b). As U x<j 
is the noncompact gauge field, the density of instanton 
pu is vanishingly small. On the other hand for V x ,j, Pv 
exhibits a hysteresis loop just like the internal energy E 
at the critical point. Vanishing of pv for c v > 0.4(0.6) 
means that the observed phase transition is the normal 
to Higgs-SC phase transition. As in the compact U(l) 
gauge case [9] , adding small but finite positive G? m -terms 
stabilizes the sign of (UUVV) as (UUVV) > in the 
Higgs-SC phase. This SC phase corresponds to the ex- 
tended s-wave, because on-site amplitude of the Cooper 
pair is zero whereas expectation values of V x .j on links 
( x ij) (j = 1,2,3) have the same sign under the gauge- 
fixing condition ip x — in Eq.(9) [13]. 




FIG. 3. Results for ^4gl with c m = 1, d m = 0, c u = 0.5 
(L = 16). (a) Internal energy exhibiting hysteresis loop, (b) 
Densities of instanton. pu (pv) is the instanton density of the 
noncompact gauge field Ax ,j (compact gauge field V x j). 
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FIG. 4. Phase diagram of AaL with c m = 1, d m = 0. Lo- 
cations of the phase transition points are determined by those 
of center of the hysteresis loops. "Normal" phase denotes the 
Coulomb (confinement) phase of U x ,j (V x j). 
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FIG. 5. Results for Agl with c m = 1, d m = 0, c„ = 0.5. 

(a) Gauge-boson mass (L = 16) measured for increasing c v . It 
exhibits a sharp discontinuity at c v = 0.6, the phase transition 
point, (b) Gauge-boson mass for decreasing c v . 

In Fig. 4, we show the phase diagram obtained from 
the measurement of the internal energy. In Fig. 5(a) and 

(b) , calculations of the mass of the gauge boson A. x j 
are given. Contrary to the AHM in Fig.l, the gauge- 
boson mass exhibits a sharp discontinuity and acquires 
nonvanishing value at c v = 0.60(0.45). 

We also studied phase transition with c u , c v fixed and 
c m varied, and found that a first-order phase transition 
occurs at a critical coupling of c m . Measurement of E 
is given in Fig. 6 for the c u = c v = 1 case. Density of 
instantons exhibits similar behavior. 
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FIG. 6. E as a function of c m for c u — c v = 1. 



Let us turn to the d m < case. In the large-c u limit, 



U x ,j ~ e %VxJti e lVx and the <i m -terms in Aql prefer con- 
figurations like (e-'^+'y^i^e'^w) < (i ^ j). Then 
under the gauge-fixing condition ip x = 0, it is expected 
that the expectation value of the Cooper-pair field (V x ,i) 
changes its sign under a 7r/2 rotation in a plane. Though 
some of the d-wave SC materials have a layered structure, 
we first consider the 3D isotropic case and set d m = —0.8. 

Phase structure was studied by means of the MC sim- 
ulations as before and found that there exist phase tran- 
sition lines. The internal energy E shows hysteresis loop 
at critical points as in the previous case. We show the 
calculations of E for certain places in the c„ — c v plane 
in Fig. 7. Phase diagram obtained by the measurement 
of E and C is given in Fig. 8. The gauge-boson mass Mq 
exhibits a discontinuity at phase transition points as in 
the previous case. 
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due to the frustration, is not observed in this case and 
the first-order phase transition to the Higgs-SC phase is 
observed instead. We also verified that critical line of the 
normal-SC phase transition exists as in the isotropic case. 
However, discontinuity in the gauge-boson mass Mq at 
critical points becomes smaller than that in the isotropic 
case. 

In conclusion, we studied the noncompact U(l) lattice 
gauge model with link Higgs field that is a GL theory for 
the unconventional SC including the extended-s, d-wave 
and also ferromagnetic p-wave SC's. By means of the MC 
simulations, we clarified the phase structure. There exist 
first-order phase transitions from the normal to Higgs-SC 
phases. We also observed that as the anisotropy of the 
layered structure is getting larger, signal of the first-order 
phase transition is getting weaker. 
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FIG. 8. Phase diagram for c m = 0, d m = —0.8. Loca- 
tions of the phase transition points are determined by those 
of center of the hysteresis loops. 



Besides the normal and Higgs-SC phases, we found 
that there exists an exotic phase that we call staggered 
state (SS). Existence of a similar phase has been ob- 
served for the compact gauge model [9]. It stems from 
the fact that in 3D there are no configurations that sat- 
isfy V x ,iV2 j < (i 7^ j) for all i and j simultaneously. In 
other words, the <i m -terms cause frustrations. E and the 
instanton density exhibit the first-order phase transition 
at the boundary of the normal and SS states. However, 
we also found that as the system size is getting larger, 
signal of the phase transition at the boundary is getting 
weaker. It is possible that finite region of the SS disap- 
pears for the infinite system-size limit. 
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FIG. 9. Results for Aql with c m = 0, d m = -0.8(-0.4) 
for intralayer(interlayer) coupling and c u — 0.5 (L — 24). (a) 
E as a function of c v for c u — 0.5. (b) Measurement of Mq 
(average of Mg in three directions) as a function of c v for 
increasing c„. Values of Mg are much smaller than those in 
the 3D isotropic case. 

Finally let us turn to anistropic cases and study how 
the phase structure changes due to the layered structure. 
To this end, we put different values for the interlayer and 
intralayer d m 's. Numerical results for d m = — 0.8(— 0.4) 
for the intralayer (interlayer) coupling with c u = 0.5 are 
given in Fig. 9. The SS, which exists in the isotropic case 
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